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We give a theoretical study of unusual resistive (dynamic) localized states in anisotropic Josephson 
junction ladders, driven by a DC current at one edge. These states comprise nonlinearly coupled 
rotating Josephson phases in adjacent cells, and with increasing current they are found to expand into 
neighboring cells by a sequence of sudden jumps. We argue that the jumps arise from instabilities in 
the ladder's superconducting part, and our analytic expressions for the peculiar voltage (rotational 
frequency) ratios and I-V curves are in very good agreement with direct numerical simulations. 

PACS numbers: 05.45.-a, 74.50. +r 



I. INTRODUCTION 

The subject of large-amplitude anharmonic dynamics 
in lattices has received widespread attention over the past 
decade. In particular, intense theoretical focus has cen- 
tered on so-called intrinsic localized modes (ILMs), also 
known as discrete breathers, with the result that many of 
their properties are now well understood These ex- 
citations result from the interplay between nonlinearity 
and discreteness, and they can be highly localized in per- 
fect lattices, with or without external driving. They can 
occur in a variety of different lattices: recent experiments 
have reported vibrational ILMs in a quasi-lD charge- 
density wave system spin- wave ILMs in a quasi-lD 
antiferromagnetic system and discrete breathers in 
Josephson junction (J J) ladders |^,||. 

The latter systems are noteworthy, in that arrays of 
coupled JJs have served for many years as reliable labo- 
ratory systems for studying diverse nonlinear phenomena 
1^. The nonlinear dynamics are particularly rich. A sin- 
gle "small" JJ subject to an applied constant DC bias 
current can be mapped onto the problem of a damped 
pendulum driven by a constant torque, with the dy- 
namical degree of freedom being the Josephson phase 
difference (Q. There are thus two qualitatively differ- 
ent states, namely a static (superconducting) state and 
a dynamic (whirling or resistive) state, with the latter 
producing a readily measured voltage V (x if across the 
whirling junction. When several junctions are assembled 
to form a regular array, such as the ladder shown in Fig. |^, 
they become inductively coupled. In the coupled system, 
junctions in the superconducting state can also exhibit 
steady-state librations, when JJs in the whirling state 
are present. In view of the mapping onto the pendulum 
problem, J J ladders share features with lattices of nonlin- 
early coupled electric dipole rotors, driven by an external 
monochromatic AC electric field |^ , but with the impor- 
tant simplification that they can be driven with purely 
DC bias currents. 

Figure |l| sketches an anisotropic JJ ladder consisting 
of small JJs of two types, "horizontal" and "vertical," 



which are, respectively, perpendicular and parallel to the 
applied bias current (arrows) . The anisotropy arises from 
the different areas of the horizontal and vertical junctions 
and is characterized by the parameter rj = A^/A^ — 
Ic/Iq, where and are the the critical currents for 
each type of junction. 
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FIG. 1. Sketch of a ID Josephson junction ladder. Solid 
arrows: applied DC edge current. Dashed arrows: bulk cur- 
rent. 



References yj5|,g| report experimental observations of 
various discrete breathers in ladders driven by a homo- 
geneously applied DC bias current, represented by the 
dashed arrows in Fig. For these states, the local- 
ized voltage patterns have a simple structure involving 
only two nonzero steady-state voltages (rotational fre- 
quencies). The breathers were found to be stable in the 
limit of small coupling {rj ^ 1) and for bias currents 
-^oxt ^ Ic- For the case of large 2-D J J arrays subject 
to a homogeneously applied DC bias, more complicated 
inhomogeneous states, with meandering voltage patterns, 
have also been reported |10|. 

Here, we study the dynamics of a JJ ladder with an 
external DC bias current applied at only one edge (solid 
arrows. Fig. |^). For increasing bias {Icxt ^ Ic) and over 
a wide range of anisotropics, we find by direct numeri- 
cal simulations that the dynamic state expands into the 
ladder one cell at a time, by a sequence of abrupt jumps. 
This behavior is in marked contrast to the well-known 
cases of long JJs and J J parallel arrays (ji = oo), where 
the entire system abruptly switches to the resistive state 
at a particular value of the DC bias. It is also different 
than the breather case, since all of the junctions within 
the boundary of this localized dynamic state whirl, and at 
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each expansion the number of different frequencies (volt- 
ages) grows. The sequence of I-V characteristics and 
threshold currents can be modeled analytically, yielding 
very good agreement with the numerical results. 



II. NUMERICAL SIMULATIONS 

We consider a ladder with a large but finite number 
of cells N. The ladder's state is specified by the time- 
dependent Josephson phases {ipi}, and {"ipi} for the 
vertical, upper horizontal and lower horizontal junctions, 
respectively, where i denotes the cell. We have found in 
our simulations that the symmetry condition ipi = —ipi 
holds for the phenomena to be discussed here. The ladder 
dynamics are then determined by the coupled nonlinear 
equations of motion obtained in Refs. ||ri| and |p^ : 

PL 

i = 2,...,iV-l, (1) 

i(V'i) = — j-Cv'i - <y5j+i - 27/-,), i = l,...,N, 

VPl 

where the operator L{ip) = ip + aip + sm{(p). The equa- 
tions for the vertical junctions at « = 1 and N arc 



H'Pl) = 71 + ^("^2 -(/?!+ 2V'i), 

PL 
PL 



(2) 



Equations and (^) describe each junction within 
the resistively and capacitively shunted junction (RCSJ) 
model , and the unit of time is the inverse of the plasma 
frequency = \/2eIc/ Ch. Since each junction's critical 
current and capacitance scale with the area, loj is inde- 
pendent of the anisotropy parameter 77, as is the effec- 
tive damping constant a = l/{ujjRC). The normalized 
bias current 7^ is defined as li.cxt/Ic- The inductive cou- 
pling between junctions is determined by the parameter 
Pl = 27rL/"/$o, where L is the self-inductance of a sin- 
gle cell and $0 — hc/2e is the elementary flux quantum. 
Coupling beyond that described by Pl is not included. 

We performed direct numerical integration of the equa- 
tions of motion for ladders with N = 20 cells, using a 
fifth-order Gear predictor-corrector algorithm [[l3| , for a 
range of anisotropics: ry — 0.5, 1.0, 2.0, 3.0, and 5.0. The 
arrays were underdamped, with a = 0.1, and we used a 
moderate value of the coupling parameter = 0.5. The 
external DC bias was applied at one edge, i.e. 71 — 7 and 
all other 7^ = 0. To simulate the I-V curves, we started 
with all phases at zero and gradually increased the exter- 
nal bias 7 from zero to 50, in increments of 0.005. When 
junctions were present in the whirling state, the MD time 
scale was set by the time-average period of the fastest ro- 
tating phase. For a given value of gamma, we waited for 



at least 100 of these reference periods before averaging, in 
order to avoid transients, following which we computed 
the time-averages {(fi) and (tpi) over at least 100 addi- 
tional reference periods. These averages are proportional 
to the average voltages across the junctions. The current 
was then incremented, with the initial phase configura- 
tion being that from the preceding MD time step. In all 
runs, the time step was 1/200 of the reference period. 
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FIG. 2. Numerically simulated I-Vi curves for the first four 
cells in a 20-celI JJ ladder with anisotropy r\ = 2. The top and 
bottom panels give the time-average frequencies (dimension- 
less voltages) of the Josephson phases for vertical and hori- 
zontal junctions, respectively. The addition of each new cell 
to the dynamic state is accompanied by abrupt jumps for the 
other whirling phases. The dimensionless DC bias threshold 
currents for the jumps are 7 = 2.040, 4.075, 10.20, and 26.52. 
Extending the current to 50 yielded no additional cells. 

Our simulated 1-Vi curves for an anisotropy oi r] — 2 
are shown in Fig. ^ The most striking finding is the oc- 
currence of extremely sharp voltage jumps. At each of the 
corresponding threshold currents 7^'"', a new cell is added 
to the ladder's dynamic state. With the applied current 
below the first threshold 7^'"", all junctions are in a static 
(zero voltage) state. When 7 exceeds 71'"^, the first verti- 
cal junction and its adjacent top and bottom horizontal 
junctions abruptly switch into the dynamic state, with 
all other junctions remaining in the zero voltage state — 
the rotating phases are confined to the first cell. As the 
bias is increased further, all three average voltages for 
this 1-cell dynamic state increase linearly until the next 



2 



threshold current 72'"' is reached, at which point the dy- 
namic state suddenly expands into the second cell, ac- 
companied by sharp changes of the voltages in the first 
cell. This process continues, yielding successive transi- 
tions from n-cell dynamic localized states to (n + l)-cell 
dynamic states. The distribution of threshold currents 
and voltage ratios depends on the ladder's anisotropy. 
Over the range < 7 < 50, the ladder reached a 3-cell 
state for 7/ = 0.5 and 1, a 4-cell state for 77 = 2, and a 
5-cell state for 77 = 3 and 5. In the following, these states 
will be termed n-cell edge states. 
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FIG. 3. Instantaneous images of MD-simulated Josephson 
phase distributions in the left half of the rj = 2 ladder underly- 
ing Fig. ^ for different applied edge currents. The dashed cir- 
cles denote junctions in the whirling state. Panel (a): phases 
immediately before the appearance of the 1-cell edge state. 
Panels (b) through (e): phases just after expansion of the 
edge state into 1, 2, 3, and 4 cells, respectively. Note the 
vortex state in the superconducting part of the ladder. 

An ri-cell edge state is in striking contrast to an n-cell 
breather. The breather occurs away from the ladder's 
edge and is homogeneously driven by same DC bias cur- 
rent (/oxt ^ Ic) applied at every cell, whereas the edge 
states are driven by a DC bias (/cxt ^ Ic) applied at 
just one edge. The edge states have a richer internal 
structure — all of the junctions within an edge state are in 
a nonzero voltage state (see Fig. ||), whereas in a breather 
state, all of the horizontal junctions are in the zero volt- 
age state except for those on the breather's boundary 
Moreover, all of a breather's vertical junction 
phases rotate at the same average frequency, whereas 



the ?i-cell edge state exhibits a peculiar distribution of 
average frequencies. This frequency (voltage) distribu- 
tion depends on both n and the ladder's anisotropy. For 
example, our simulations for the rj — 2, 3-cell edge state 
of Fig. |(d) yield the ratios given in second column of 
Table |. Comparison with the third column shows that 
they are in very good agreement with analytic predictions 
derived below. 

TABLE I. Average frequency (voltage) ratios for 3-cell 
edge states. The MD ratios are for the rj — 2, 3-cell edge 
state of Fig. ^(d), and the predicted ratios were calculated 
from Eqs. (|) and (|). 



Ratio 


MD 


Predicted 


V / h 


2.667 
8.006 
2.499 
5.017 
2.005 


(3J7^-f 8r?-H4)/2r;(l-f r?) =8/3 
(377^ +8ri + i)/rf = 8 
(77" +6»7 + 4)/»7(2 + »?) = 5/2 
{rf +Qrj + 4)/if = 5 
2 



The superconducting state forming ahead of the n-cell 
edge state is also unusual. Fig. || gives snapshot images 
of the Josephson phase distribution for several values of 
the applied DC current bias, for the anisotropy rj — 2. 
In panel (a), the current is just below the first threshold, 
and one sees a single Josephson vortex in the supercon- 
ducting part of the ladder. The remaining panels (b)-(e) 
show the phases just after a new cell is added to the dy- 
namic state. At the threshold currents 7^^"", the super- 
conducting state becomes unstable, and the vortex jumps 
into the next cell as the edge state expands. Our simu- 
lations show that in general the superconducting state is 
sensitive to the anisotropy. Thus for rather small values 
of 77 ^ 1, there are no vortices trapped in the super- 
conducting part of the ladder over the range of currents 
studied. For these cases, the Josephson phases of the 
vertical junctions in the superconducting portion of the 
ladder simply decrease with distance from the boundary 
of the resistive portion, corresponding to the Meissner 
state of the superconductor. With increasing anisotropy, 
single Josephson vortices appear in the superconducting 
portion, as in Fig. ^. For large anisotropies (77 ~ 5) 
more complex vortex trains are observed, and we also 
find that the penetration of the edge state changes the 
nature of the superconducting vortex state, rather than 
simply pushing it ahead as for rj = 2. A detailed discus- 
sion of the superconducting state will be given elsewhere 
(Ref. 0). 



III. THEORETICAL ANALYSIS 

The unusual voltage distributions in the n-cell edge 
states can be explained analytically by making use of 
Kirchhoff's laws, applied to the time-average currents 
(normalized to /") and corresponding dimensionless volt- 
ages in each cell. The key assumption is the coexistence 
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of the resistive and superconducting states in different 
portions of the ladder. For a cell i within an n-cell edge 
state, current conservation gives + /f = //Li, while 
the voltage condition is /" — — T"_^^ = 0. Combining 
these yields an equation for just the horizontal currents: 



100 



rh 



2 1 



0. 



(3) 



This equation and the two from which it was derived 
apply to all cells i in 1 < i < n, provided we define 
Iq = 7, r^^-^ = 0, and /^^^ = in order to take 
proper account of the n-cell edge state's boundaries. 
Equation (||) is readily solved by substituting = A' , 

which yields two roots, namely A = 1 + ;^ + ^/(l + — 1 

and 1/A. Hence the general solution of Eq. (^) is = 
CiA' + C2A~', where the constants Ci and C2 are obtained 
from the above definitions of Iq and I^+i- With the 
horizontal currents thus determined, the vertical currents 
may be computed from If = I^_-^ — I^. The currents are 
then converted into the average voltages via Vj" — 1^ /a 
and Vl^ = I^/iai]). The resulting voltage distribution 
within an n-cell edge state is (1 < z < n): 



7(l-A)(A'-i -A^"+i-') 
a(A2«+i + 1) ' 



7(A* + a2"+i-*) 
ari{X 



2n+l 



1) 



(4) 



(5) 



Equations (||) and (||) give the predicted voltage ratios 
in Table |, which for rj = 2 are seen to be in very good 
agreement with our MD simulations. Indeed, we find 
that for all of the values of rj studied, the predicted I-V 
curves are in very good agreement with the MD curves, 
such as those of Fig. ^. Only the values of the current 
thresholds for the jumps are left undetermined by these 
equations. 

We can also predict the distribution {7^'""} of threshold 
currents for each i] by assuming that the superconducting 
state associated with the (n — l)-cell edge state becomes 
unstable and converts to the n-cell edge state when the 
current I^ reaches a depinning current I^p, which we take 
to be independent of n. This yields an expression for the 
threshold currents 



.thr 



I 



dp 



cosh[(n — ^)lnA] 
cosh (| In A) 



(6) 



The ratios if'' /I dp predicted by Eq. (|) for n 2, 3, 4 
and 5 are plotted versus rj in Fig. ^. To compare with the 
MD results, we fit I dp to the first observed MD threshold 
current for each r?, namely 7^^"' = 1.295, 1.510, 2.040, 
2.525, and 4.010, for ij = 0.5, 1.0, 2.0, 3.0, and 5.0, re- 
spectively. The resulting MD ratios are shown by the 
circles in Fig. and agree well with the predictions. 
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FIG. 4. Threshold current ratios 7^'"^//dp. Curves: pre- 
dicted ratios vs. anisotropy 77, from Eq. (^. Starting at the 
bottom, the curves are for n = 2, 3, 4 and 5. Circles: numer- 
ical ratios from the MD simulations. 



IV. HYSTERESIS 

Despite their rich structure of frequency ratios, the n- 
cell edge states are found to be highly stable in our sim- 
ulations, for the case of increasing current. However, we 
also find notable hysteresis effects when the simulations 
are started with an n-cell edge state for large n and the 
applied DC current is gradually decreased to zero. Fig- 
ure ^ is representative of the hysteretic behavior encoun- 
tered. The threshold currents and stability properties for 
the sequence of down-conversions {n — > n — 1} are very 
different than for the increasing-current case. In partic- 
ular, we observe resonant steps, switching processes, and 
nonlinear regions in the I-V curves. We believe that all 
of these features arise from the resonant interaction be- 
tween the n-cell edge states and other excitations, both 
localized and delocalized, as will be discussed elsewhere 
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V. SUMMARY 



In summary, our numerical simulations have revealed 
unusual localized dynamic states in anisotropic JJ lad- 
ders subject to a DC bias current at one edge. Increasing 
the bias causes these states to expand by adding single 
cells in a sequence of sudden jumps, giving rise to a di- 
verse set of voltage distributions and sharp changes in 
the I-V curves. This behavior occurs for a wide range 
of parameters and should be observable through the I-V 
characteristics or by direct visualization using low tem- 
perature scanning laser microscopy techniques Pj9|JTo|] . 
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FIG. 5. Hysteresis. Shown arc the nurnericaUy simulated 
I-V curves for the vertical junctions of the first few cells in 
a 20-cell JJ ladder with anisotropy = 3. The bias cur- 
rent was first increased (dashed curves) from zero to 50 in 
increments of 0.005, resulting in a 5-cell edge state, at which 
point the current was decreased (solid curves) back to zero 
with the same increment. At each value of the current, we 
waited 400 reference periods and then averaged the voltages 
over an additional 200 periods. Only a small current range 
is shown, in order to illustrate clearly the differences between 
the increasing- and decreasing-current cases. The latter case 
has a different sequences of threshold currents and exhibits 
nonlinear regions and small steps not associated with the dis- 
appearance of the edge states. 
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